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Abstract 

Carbonaceous chondrites are usually 
believed to be the primary constituents of 
near-Earth asteroids and Phobos and 
Diemos, and are potential resources of 
fuels which may be exploited for future 
planetary missions. In this work we have 
calculated the nature and abundances of 
the major volatile bearing and other 
phases, including the vapor phase, that 
should form in Cl and C2 type carbonaceous 
chondrites as functions of pressure and 
temperature. The results suggest that 
talc, antigorite ± magnesite are the major 
volatile bearing phases and are stable 
below 400°C at 1 bar in these chondritic 
compositions. Simulated heating of a 
kilogram of C2 chondrite at fixed bulk 
composition between 400 and 800°c at 1 bar 
yields about 135 gm of volatile, which is 
made primarily of H 2 o, CH 4 , C0 2 and 
CO. The relative abundances of these 
volatile species change as functions of 
temperature, and on a molar basis, H 2 
becomes the most dominant species above 
550°C. In contrast, Cl chondrites yield 
about 306 gm of volatile under the same 
condition, which consist almost completely 
of 60 wt% H 2 0 and 40 wt% C0 2 . Preliminary 
kinetic considerations suggest that equi- 
librium dehydration of hydrous phyllosili- 
cates should be attainable within a few 
hours at 600°c. These results provide the 
framework for further analyses of the 
volatile and economic resource potentials 
of carbonaceous chondrites. 

Introduction 

The near-Earth asteroids are poten- 
tial resources of volatiles which can be 
used as propellants and life supporting 
purposes for planetary missions. As 
emphasized by Lewis and Lewis (1), these 
asteroids are at times the nearest bodies 
to Earth, and many of them can pass 
between Earth and Moon. These authors 
have also shown that within any arbitrary 
three-year period, there are roughly 300 
and 90,000 launch opportunities to respec- 
tively kilometer and 100-meter sized 
asteroidal bodies as compared to 39 lunar 
launch windows. These facts and the 
potential volatile resources make the 
asteroids important candidates for 
detailed analyses for supporting future 
planetary missions. 

Cost-effective or energy-efficient 
extraction of volatiles from the asteroids 
require knowledge of the nature and modal 


abundances of the constituent minerals in 
which the volatile components are struc- 
turally bound. There has been several 
studies on reflectance spectroscopy of 
near-Earth asteroids to characterize the 
nature of the volatile bearing minerals 
(2,3). While these studies report absorp- 
tion band near 3 urn suggesting presence of 
hydrous phy llosil icates , more precise 
characterizations were equivocal. An 
obvious alternative approach to the reso- 
lution of the problem is to examine 
directly the asteroidal materials. 

The spectral and density character- 
istics strongly suggest that a significant 
fraction of the near-Earth asteroids are 
made of carbonaceous chondrites (1,4). 
The most volatile rich of these meteorites 
are what are commonly known as Cl or (Cl) 
and C2 (or CM) classes. However, the 
grain size of the volatile bearigg phases 
are often too small (100 to 1000 A) (5) to 
permit identification under optical micro- 
scope, and are best observed under high 
resolution transmission electron micro- 
scopy ( HRTEM ) . Zolensky and McSween (5) 
have recently presented a summary of vola- 
tile bearing phases identified in the 
carbonacous chondrites. 

The primary objective of this work is 
to carry out theoretical calculations to 
predict modal abundances and compositions 
of the major mineral phases, along with 
the abundance and composition of the co- 
existing vapor phase, that could develop 
in the bulk compositions of Cl and C2 
chondrites as functions of pressure and 
temperature. The results would provide 
the framework for engineering designs for 
the extraction of volatile components from 
asteroids as well as Phobos and Deimos, 
the two small natural satellites of Mars, 
which are also likely to be made of car- 
bonaceous chondrites (1). 

Theoretical Method 

Principles 

According to Duhem’s theorem in 
classical thermodynamics (6,7), the equi- 
librium state of a closed system (i.e., a 
system of fixed composition and mass) is 
completely determined if any two variables 
are fixed regardless of whether these are 
intensive, extensive, or a combination of 
both. Ganguly and Saxena (7) have recent- 
ly reviewed the various methods by which 
one may carry out the actual computation 
of the equilibrium assemblages. The 
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Table 1 Bulk chemical 
compositions of Cl and C2 
carbonaceous chondrites 




Cl 

C2 


Si 

10.40 

12.96 

Ti 

0.04 

0.06 

Al 

0.84 

1.17 

Cr 

0.23 

0.29 

Fe 

18.67 

21.56 

Mn 

0. 17 

0.16 

Mg 

9.60 

11.72 

Ca 

1.01 

1.32 

Na 

0.55 

0.42 

K 

0.05 

0.06 

P 

0.14 

0.13 

Ni 

1.03 

1.25 

Co 

0.05 

0.06 

S 

5.92 

3.38 

H 

2.08 

1.42 

C 

3.61 

2.30 

0 

45.61 

41.74 


method used in this study obtains the com- 
positions and abundances of the equilib- 
rium phase assemblages by minimizing the 
Gibbs Free Energy (G) of the system at 
fixed P-T conditions. The minimization is 
constrained to conserve the various ele- 
mental masses of the system through the 
method of Lagrangian multipliers (7). 
Because of the extremely small grain size 
of the phyllosilicates observed in carbo- 
naceous chondrites, one should consider 
the effects of surface free energy on the 
formation of these phases. As a first 
approximation, we have, however, ignored 
this effect. Therefore, the actual equi- 
librium temperatures of volatile bearing 
phases could be lower than those predicted 
by our calculations. 

System, Phases and Data Base 

Table 1 shows the bulk compositions 
of Cl and C2 carbonaceous chondrites. 
These compositions are based on analyses 
of a number of fragments of Orgueil and 


Table 2 Phases in the system 
Mg-Fe-S i-C-H-O-S considered in this work 

The sources of thermochemical data are shown by reference 
numbers within parentheses. The Heat of Formation of the Fe-end 
members of the hydrous phyllosilicates are derived in this work. 


Anhydrous Silicates and Oxides 


Olivine 

( Mg , Fe ) 2 Si0 4 

(8) 

Orthopyroxene 

(Mg,Fe)Si0 3 

(8) 

Per iclase 

( Mg , Fe )0 

(23) 

Ouartz 

sio 2 

(17) 

Cristobal ite 

S i0 2 

(17) 

Hematite 

Fe 2 ° 3 

(23) 

Magnetite 

Fe 3 ° 4 

(23) 

Hydrous and Carbonate 

Phases 


Anthophyll i te 

(Mg,Fe) 7 Si 8 0 22 (0H) 2 

(8,27) 

Talc 

(Mg,Fe) 3 Si 4 O 10 (OH) 

(8) 

Ant igor i te 

(Mg,Fe) 48 Si 34 0 85 (OH) 62 

(8) 

Chrysot i le 

(Mg,Fe) 3 Si 2 0 5 {0H) 4 

(8) 

Bruc i te 

( Mg , Fe ) ( OH ) 2 

(17,24) 

Magnesite 

(Mg,Fe)C0 3 

( 17,25) 

Sulfides and Elements 



Troil i te 

FeS 

( 17) 

Pyr i te 

FeS 2 

(17) 

Iron 

Fe 

(26) 

Sulfur 

S 

(17) 

Graphite 

C 

(17) 

Vapor Phase 

C-O-H-S 

(15,16,17) 


Species: H 2 0, C0 2 , CH 4 , CO, 0 2 , H 2 S, S 2 , S0 2 , COS 
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Murchison meteorites, respectively, and 
are summarized in Dodd (22). It can be 
easily seen from this table that the sub- 
system Mg-Fe-Si-C-H-O-S (MFCHOS) consti- 
tutes almost 96% by weight of the bulk 
compositions of Cl and C2 chondrites. In 
order to somewhat simplify the computa- 
tional problem, we have confined the G 
minimization calculation to this sub- 
system, but evidently the phases whose 
compositions form outside of this sub- 
system cannot be major constituents of 
these chondrites. 

Table 2 shows the list of phases that 
we have considered as possible crystal- 
lizing phases within the above subsystem 
along with the sources of the thermochem- 
ical data. The choice of these phases 
have been guided by the published reports 
of the mineralogy of carbonaceous chon- 
drites whose compositions lie within the 
system MFCHOS (e.g., Zolensky and McSween 
(5)). The selection of thermochemical 
data in the system Mg-Si-O-H is based on 
the critical analysis of Chatterjee (8,11) 
for consistency with both calorimetric and 
experimentally determined phase equilib- 
rium measurements. These data are similar 
to those of Berman (9). However, unlike 
those in Berman, the data are also consis- 
tent with high P-T experimental data 
involving distribution of Fe and Mg among 
coexisting silicates. 

The heat capacities (C D ) and entro- 
pies (S°) of the phases, if^not available 
in any of the recent systematizations of 
self-consistent data set, are estimated by 
stoichiometric summation of the properties 
of oxides and structurally analogous com- 
pounds, as discussed in a number of publi- 
cations (13,14). 

There are no thermochemical data for 
the Fe-end members of the phyllosi 1 icates 
listed in Table 2. The enthalpies of 
these phases are estimated according to 
the known relative properties of Fe- and 
Mg-end members of biotite, which is also a 
phyllos il icate . The data derived by 
Chatterjee ( 8 ) and Robie Hemingway 
(10) for the Mg- and Fe Z -end members 
(phlogopite and annite, respectively) 
yield an enthalpy difference of -359.566 
KJ per mole of divalent cation at 1 bar, 
298 K. We have used these values to esti- 
mate the enthalpies of formation from 
elements of Fe-end members of phyllosili- 
cates at 1 bar, 298 K from the available 
data for the Mg-end members. The results 
are as follows ( KJ/mol of Fe ) . Fe-talc: 
-1606.80; Fe-an t igor i te : -1127.53 and Fe- 
chrysotile: -1094.87. The corresponding 
values of the free energy of formation 
from elements (aG^) for Fe-talc and Fe- 
antigorite at 1 bar, 2^ K are -4477.18 
and -4784.0 KJ/mol of Fe^ , respectively. 

The validity. of the estimated thermo- 
chemical properties of the Fe-end members 
of the phyllosilicates may be tested as 
follows. The Fe-end members of talc and 


antigorite are unstable. Instead, the 
stable phases are closely analogous com- 
pounds, namely, minnesotite, which has a 

stoichiometry of Fe 27 s *36°8 6 ( 0H ) 26 as 
compared to that of Fe 27 Si 3^090 ( OH ) jo for 
Fe-talc, and greenalite with a stoicniom- 
etry of Fe4gSi3 2 Ogo ( OH ) 64 compared to 
Fe4 8 Si3 4 085 ( OH) for Fe-an t igor i te . 
Anovitz et al. ?12) have recently deter- 
mined AGf of minnesotite and greenalite 
from phase equilibrium data. Their 
results, AGf (minnesotite) = -4577.78 KJ 
and AGf (greenalite) = -4799.50 K J , 
differ in the right direction from the 
corresponding values of Fe-talc and Fe- 
antigorite, considering the relative 
stabilities of the phases. 

The stereochemical environments of 
divalent cations in biotite are quite 
similar to those in talc, both being 2:1 
(TOT) layer silicates, and partly similar 
to those in antigorite, which is an 1:1 
(TO) layer silicate (28). Thus, these 
phyllosilicates should have similar 
differences between the enthalpies of 
formation of the Fe- and Mg-end members. 
Further, since on the basis of observa- 
tional data, the predicted compositions of 
the mineral phases in carbonaceous chon- 
drites are expected to be Mg-rich, the 
effects of errors in the Fe-end member 
thermochemical properties and solid solu- 
tion model will be relatively small. This 
can be easily understood by considering 
the form of G-X curve of a stable solution 
and evaluating the effects of errors in 
above properties near Mg-terminal segment. 

The ±£ermodynamic mixing properties 
of the Fe^ - and Mg-end member components 
of the anhydrous silicates are taken from 
the self-consistent summary of Chatterjee 
(8), while those of the phyllosilicates 
are assumed to be ideal for the lack of 
any data. The latter assumption is likely 
to be approximately valid at T 600 C, 
since at these conditions biotite solid 
solution is found to behave approximately 
ideally (7). At lower temperature, the 
phyllosilicates may deviate from ideal 
solution behavior, but perhaps not too 
dissimilarly to significantly affect their 
relative stabilities. Magnesite is also 
assumed to be an ideal solution of Fe-Mg 
components. As emphasized by Ganguly and 
Saxena (7), the computed dehydration con- 
ditions are not very sensitive to errors 
in activity-composition relations since 
the devolatilization equilibria are char- 
acterized by large enthalpy changes. The 
iron oxides and sulfides are considered to 
be stoichiometric phases since the solid 
solutions of Mg are around a few percent 
at the temperatures of interest. 

The fluid properties are taken from 
Saxena and Fei (15,16) and Robie et al. 
(17). The "corresponding state method" of 
estimating P-V-T relation, as discussed 
by these authors, is extended to include 
the sulfur species. For mixing data, we 
have assumed that H 2 S behave as H 2 0, S 2 as 
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Fig. 1 Calculated modal abundances of minerals per kilogram of 
rock in (a) Cl and (b) C2 carbonaceous chondritic bulk compositions 
as functions of temperature at 2 Kb. The steps in the calculations 
are shown by symbols. 


C> 2 , and SO 2 and COS as CC^* and used the 
equivalent binary interaction parameters 
from Saxena and Fei (15). 

Results 

At 1 bar, the free energy minimiza- 
tion calculations do not yiel^ any vola- 
tile bearing phase at T > 400°C. At the 
present state of the thermochemical data, 
we do not believe that our calculations 


are reliable at lower temperatures. In- 
deed, we find results at lower tempera- 
tures which do not seem to follow the 
systematics of the results obtained at 
higher temperatures. Thus, we have calcu- 
lated the equilibrium phase assemblages at 
a pressure of 2 Kb so that the volatile 
bearing phases appear at higher tempera- 
tures. This procedure at least permits 
recognition of the major volatile bearing 
phases. Some of these phases may- be re- 
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Temperature (C) 


Fig, 2 Calculated compositions of the silicate phases in Cl and 
C2 carbonaceous chondrites as functions of temperature at 2 Kb, 


placed by more stable assemblages at lower 
temperature, but should reappear during 
heating experiment. The results, which 
are illustrated in Figure 1, show that the 
primary volatile bearing phases in the Cl 
and C2 bulk compositions are talc, anti- 
gorite ± magnesite. Talc is stable to 
around 650 C at 2 Kb, and dehydrates to 
orthopyroxene (OPx) at higher temperature 
according to Tc + Ol » 5 OPx + H 2 0. Anti- 
gorite and magnesite crystallize at T < 
550°c, but magnesite is absent in C2 bulk 
composition at T >_ 400°C. FeS is a stable 
phase at all temperatures investigated in 
this work. We, however, did not consider 
the formation of sulfate or pyrrhotite. 
Magnetite becomes a stable phase at 
T <_ 500°C. Elemental sulfur found in Cl 
chondrite is invariably associated with 
pyrrhotite (Fei_ x S) (5). This suggests 
decomposition of troilite (FeS) to pyrrho- 
tite and sulfur. 

Figure 2 shows the compositions of 
the f er romagnes ian phases as function of 
temperature at 2 Kb for the Cl and C2 bulk 
compositions. Calculations at 1 bar also 
yield similar compositions. The observed 
compositions of these phases in the car- 
bonaceous chondrites vary widely in Fe/Mg 
ratio. One way to achieve this diverse 
Fe/Mg ratio is by aqueous alteration which 
seems to be a cha racter is Ly: of the car- 
bonaceous chondrites. Fe^ fractionates 
very strongly into an aqueous phase rela- 
tive to the silicates (18). Thus, varying 
degrees of equilibration with an aqueous 
phase could lead to a spectrum of Fe/Mg 
compositions of the ferromagnesian sili- 


cates. Aqueous alteration may not be the 
sole reason, but is likely to be an impor- 
tant contributing factor in the develop- 
ment of variable Fe/Mg ratio of the ferro- 
magnesian silicates in carbonaceous chon- 
drites. 

The abundances of volatile species in 
equilibrium with the solids per kilogram 
of total mass of C2 chondrites are illus- 
trated in Figure 3. Figures 3a and 3b 
illustrate the molar and mass abundances, 
respectively, as functions of temperature 
at 1 bar. It should be noted that the 
total mass of the volatiles do not change 
above 400 U C, which means that under equi- 
librium condition C2 chondritic material 
will completely devolatize if heated above 
400°c. The total yield of the volatile is 
about 14% of the initial mass of C2 chon- 
drites, The change in th^ abundance of 
various species above 400 C is due to 
homogeneous reactions within the fluid 
phase. Above 550 C, hydrogen has the 
highest molar concentration in the vapor 
phase. Hydrogen is a very good propellant 
and reducing agent which may thus be pref- 
erentially extracted by heating the C2 
chondritic material to T > 550 c. 

In contrast to the C2 chondrites, the 
volatiles given off by Cl material above 
400° C consist essentially of H 2 0 and C0 2 
at 1 bar pressure. Per kilogram of mate- 
rial, the total yield of volatile is about 
306 grams of which nearly 60 wt% is H 2 0. 

The above results are subject to 
considerable uncertainties owing to the 
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Fig. 3 Calculated equilibrium (a) molar and (b) mass abundances of 
the volatile species per kilogram of C2 chondritic material. 


approximate nature of the thermochemical 
properties of the phy 1 los il icate solid 
solutions, and also since we have ignored 
the effects of surface free energy on the 
formation of extremely fine grained crys- 
tals. Further, carbonaceous chondrites of 
any given type do not have homogeneous 
compositions, but instead probably repre- 
sent agglomeration of materials formed at 
different conditions at different parts of 
the solar system. However, these calcula- 
tions are still useful in providing an 
idea of the relative stabilities, as well 


as upper stabilities of the volatile bear- 
ing phases and their average abundance. 
It is interesting to note that recent 
reflectance spectroscopic studies of car- 
bonaceous chondrites seem to corroborate 
the predictions made above as to the 
nature of the most abundant hydrous phases 
(Gaffey, pers. comm.). 

Kinetics of Devolatilization 

The equilibrium calculations pre- 
sented above provide a framework for the 


Pressure (kbar) 
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7TaIc - 3Anth + 4Qz + 4H20 



Fig, 4 Dehydration of talc as functions of time, temper- 
ature and extent of overstepping above the equilibrium 
dehydration boundary at 1 Kb. Volume fraction is that of 
talc remaining after reaction. Experimental data points, 
shown by asterisks, are from Greenwood (19). 



Fig. 5 Equilibrium relations of talc, 
antigorite and brucite along with other 
minerals in the system MgO-S i02~^2^ * 
Reproduced from Evans and Guggenheim (21). 


evaluation of the volatile as well as some 
mineralogical resource potential of Cl and 
C2 carbonaceous chondrites. The actual 
extraction of volatiles as a function of 
temperature and heating rate, however, 
depends on the kinetics of the devolatili- 
zation process. For kinetic analysis, it 
is essential to identify the minerals in 
which the volatiles are most likely to be 
structurally bound. The above results 
suggest that talc, antigorite and magne- 
site are likely to be the major volatile 
bearing phases in carbonaceous chondrites. 

Except for pure Mg-talc, there are no 
kinetic data for the other volatile bear- 
ing phases. The dehydration kinetics of 
talc was determined by Greenwood (19). At 
a given temperature, the rate o| reaction 
depends on the rate constant, K * and the 
extent of departure from the equilibrium 
boundary. This is expressed in transition 
state theory (20) according to the rela- 
tion Rate =* K ( 1-exp ( nAG/RT ) ) , where n is 
a constant, commonly assumed to be unity, 
and AG is the Gibbs free energy of the 
reaction at the P,T condition of interest. 
Figure 4, which is constructed from Green- 
wood's experimental data, illustrates the 
dehydration kinetics of talc as functions 
of time, temperature and extent of over- 
stepping (AT) above the equilibrium dehy- 
dration boundary at 1 Kb. It is evident 
from this figure that essentially complete 
dehydration of talc can be achieve^ within 
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6 hours at about 160°C above the equilib- 
rium boundary. 


and the early solar system. University 
of Arizona Press, Tucson. 


The stability fields of the pure Mg- 
end members of various phyllosi 1 ica tes of 
interest in this study are summarized for 
the condition of P tptal = P H20 in Fil 3 ure 
5, which is reproduced from Evans and 
Guggenheim (21). According to our re- 
sults, talc and other phy llos i 1 ica tes com- 
pletely dehydrate at _>_ 400°C at 1 bar in 
Cl and C2 bulk compositions. This is com- 
patible with the phase relations shown in 
Figure 5. Thus, based on the results 
shown in Figure 4, one would anticipate 
that essentially complete dehydration of 
all phyllosil icates in Cl and C2 chondrit- 
ic material is likely to be achieved ig a 
matter of a few hours at T >_ 600 c. 
Detailed studies of the devolatilization 
kinetics of both phy 1 los ilicates and car- 
bonates as functions of temperature, fluid 
composition, and grain size are currently 
under way to evaluate the practical limi- 
tations of the extraction of volatiles 
from carbonaceous chondritic materials. 
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